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Properties of dense quark matter in and out of chemical equilibrium 
are studied within the SU(3) Nambu-Jona-Lasinio model. In addition to 



' the 4-fermion scalar and vector terms the model includes also the 6-fermion 

00 ' 

flavour mixing interaction. First we study a novel form of deconfined matter, 
' meso-matter, which is composed of equal number of quarks and antiquarks. 

o' 

It can be thought of as a strongly compressed meson gas where mesons are 
Q_i' melted into their elementary constituents, quarks and antiquarks. Strongly 

bound states in this quark-antiquark matter are predicted for all flavour com- 
binations of qq pairs. The maximum binding energy reaches up to 180 MeV 
per qq pair for mixtures with about 70% of strange (ss) pairs. Equilibrated 

5-H 



baryon-rich quark matter with various flavour compositions is also studied. 
In this case only shallow bound states appear in systems with a significant ad- 
mixture (about 40%) of strange quarks (strangelets) . Their binding energies 
are quite sensitive to the relative strengths of scalar and vector interactions. 
The common property of all these bound states is that they appear at high 
particle densities when the chiral symmetry is nearly restored. Thermal prop- 
erties of meso-matter as well as chemically equilibrated strange quark matter 
are also investigated. Possible decay modes of these bound states are dis- 
cussed. 
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I. INTRODUCTION 



The main goal of present and future experiments with ultrarelativistic heavy ions is to 
produce and study in the laboratory a new state of strongly interacting matter, the Quark- 
Gluon Plasma (QGP). In achieving this goal one is facing two major problems. First, the 
phase structure of QCD is not fully understood yet. Second, the matter evolution in the 
course of an ultrarelativistic heavy-ion collision may be out of thermodynamical equilibrium. 

Most calculations of the QCD phase diagram are made under the assumption of thermal 
and chemical equilibrium. The QCD lattice calculations at zero baryon chemical potential 
reveal a second order phase transition or a rapid crossover at temperatures around 140- 
160 MeV. Recent calculations based on different effective models |]-§J show the possibility 
of a first order phase transition at finite baryon densities and moderate temperatures. The 
predicted phase diagram in the T — \i plane contains a first order transition line terminating 
at a critical point (T c , /x c ) with T c ~ 120 MeV and a finite /i c [@||. Possible signatures of this 
point in heavy-ion collisions were discussed recently in Ref. ||. The problem is, however, 
that the matter produced in central heavy-ion collisions at very high energies has rather 
low net baryon density. If such a baryon-free matter would expand following (locally) an 
equilibrium path, it would miss the first order transition line. In this case no clear signatures 
of the phase transition would be observed. 

Moreover, it is most likely that the matter evolution in ultrarelativistic heavy-ion colli- 
sions does not follow thermodynamical equilibrium. The reason is that the matter produced 
in such collisions expands very fast. As has been already observed in heavy-ion experiments 
at the SPS energies (see e.g. ||), the expansion velocities along the beam direction are close 
to the speed of light and the transverse velocities are close to 0.5 c. A strong collective 
expansion of matter is expected also at RHIC and LHC energies. Under such conditions one 
may expect significant deviations from thermodynamical equilibrium , especially from the 
chemical equilibration. This may happen on both the partonic and the hadronic stages of 
the reaction. 
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One can mention at least two mechanisms which may lead to the deviation from chem- 
ical equilibrium on the partonic stage. First, it is believed that multiple color strings are 
produced initially in hard nucleon-nucleon collisions. Later on they decay via the Schwinger 
mechanism into quark-antiquark pairs whose abundances are determined by quark masses 
and a string tension constant. At this stage, multiplicities of secondary quarks and anti- 
quarks may be different from their values in thermodynamical equilibrium. Second, simple 
fits of QCD lattice data HU indicate that gluons may acquire a large effective mass around 
the deconfinement transition point. Therefore, even if the ideal QGP were created at some 
intermediate stage of the reaction, gluons would subsequently decay into lighter quark- 
antiquark pairs. Hence, the abundances of different quark species may deviate strongly 
from their equilibrium values. In particular, an overpopulation of the light quark-antiquark 
pairs may be expected. The phase diagram of such a chemically nonequilibrated matter may 
be very different from the predictions based on the equilibrium concepts. 

Another motivation to study chemically nonequilibrated quark-antiquark systems comes 
from the hadronic spectroscopy. It is well known that some mesonic resonances do not 
fall into the classification scheme based on the constituent quark model. They cannot be 
interpreted as conventional qq bound states. Rather they could be either bag-like (qqqq) 
states or meson-meson bound states (qq — qq). Well known examples include the /o(980) 
(KK bound state close to the threshold), the /i(1420) (KK*), the / (1500) and / 2 (1565) 



(lulu and pp) etc. [nj. A legitimate question is: what will happen if more and more qq pairs 
will be put together? Such multi-gg systems might be even more bound due to the reduced 
surface energy as compared with the bulk one. The mesonic substructure will most likely 
melt away and such states will look like multi-gg bags. We call this hypothetical state of 
matter as "meso-matter" and its finite droplets at "mesoballs". The analogous state of 



hadronic matter, bound multipion droplets, has been considered in Ref. [11 



Since the direct application of QCD at moderate temperatures and nonzero chemical 
potentials is not possible at present, more simple effective models respecting some basic 
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symmetry properties of QCD are commonly used. One of the most popular models of 
this kind, which is dealing with constituent quarks and respects chiral symmetry, is the 
Nambu-Jona-Lasinio (NJL) model fi2| , [13|| . In recent years this model has been widely used 



for describing hadron properties (see reviews [|14],[n|]), phase transitions in dense matter 
|T^n|f|,f2D| and multiparticle bound states . 



In the previous paper [^5j we have used the NJL model to study properties of the 
quark-ant iquark plasma out of chemical equilibrium. In fact, we considered a system with 
independent densities of quarks and antiquarks. To our surprise, we have found not only 
first order transitions but also deep bound states even in the baryon-free matter with equal 
densities of quarks and antiquarks. The consideration in Ref. [[25[ was limited to systems 
composed of either light (u, d) or strange (s) quarks and antiquarks. In the present paper we 
extend the model to arbitrary mixtures of light and strange quarks. The emphasis is put on 
investigating the possibility of bound states in such systems at various flavour compositions 
of quarks and antiquarks. As a special case we consider the baryon-rich quark matter in 
chemical equilibrium, in particular, the possibility of bound states in strange quark matter, 
i.e. strangelets. Thermal properties of meso-matter and strange quark matter are also 
studied. 

The paper is organized as follows: in Sect. II a generalized NJL model including flavour- 
mixing terms is formulated in the mean-field approximation. Then in Sect. Ill the model is 
used to study the bound states in qq systems with different strangeness contents. The model 
predictions for strangelets are discussed in Sect. IV. The characteristics of bound states at 
zero temperature are summarized in Sect. V. Effects of finite temperatures are considered in 
Sect. VI. Possible decay modes of new bound states are discussed in Sect. VII. Main results 
of the present paper are summarized in Sect. VIII. 
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II. FORMULATION OF THE MODEL 



We proceed from the SU(3)-flavour version of the NJL model suggested in Ref. [26|. The 
corresponding Lagrangian is written as (% = c = 1) 



j=0 

8 

-GvY, 

j=0 

-K[det f 



l-7s)^ +det / U(l + 75 ) 



(1) 



Here ip is the column vector consisting of three single-flavour spinors ?/>/, f — u,d,s, 



Ai,...,As are the SU(3) Gell-Mann matrices in flavour space, Ao = y2/3J, and 
rho = diag(mou, mod, mos) is the matrix of bare (current) quark masses. At mo = this 
Lagrangian is invariant with respect to SUl(3) <S> SUr(3) chiral transformations. The sec- 
ond and third terms in Eq. ([]]) correspond, respectively, to the scalar-pseudoscalar and 
vector-axial-vector 4-fermion interactions. The last 6-fermion interaction term breaks the 
C/a(1) symmetry and gives rise to the flavour mixing effects. In particular, this term is 
responsible for the large rj' mass ||27| . 

In the mean-field approximation the Lagrangian ([[]) is reduced to 



£ mfa = 53^7 {0- m f) 4>f 
f 



2 ' 1^11/^7 

/ 



(2) 



where Jp 



i^oGvPvf and 



Psf = < $ft>t > > 

p Vf = < Ipfjo^f > 



(3) 
(4) 



are scalar and vector densities of quarks with flavour / . Angular brackets correspond to the 
quantum-statistical averaging. The constituent quark masses, mf, are determined by the 
coupled set of gap equations 



m f = m f - G s Psf + 2K Yi Psf ■ ( 5 ) 

The NJL model is non-renormalizable, because its coupling constants have non-trivial 
dimensions: Gs, Gy oc [mass] -2 and K oc [mass] -5 . As a consequence, the contribution of 
negative energy states of the Dirac sea are divergent, and one must introduce an ultraviolet 
cut-off. In this respect the NJL model is an effective model, aimed at describing the non- 
perturbative regime of QCD at low energies. Following common practice, we introduce the 
3-momentum cut-off A to regularize divergent integrals. The structure of the fermionic 
vacuum within the NJL model is shown schematically in Fig. 1(a). Only "active" levels of 
the Dirac sea, i.e. with p < A are included in calculations. 

The model parameters mof,Gs,K,A can be fixed by reproducing the observed masses 
of 7r, K , and rf mesons as well as the pion decay constant /„.. As shown in Ref. ||26| , a 
reasonable fit is achieved with the following values: 

m 0u = m od = 5.5 MeV, m 0s = 140.7 MeV, (6) 
G s = 20.23 GeV~ 2 , K = 155.9 GeV~ 5 , A = 0.6023 GeV. (7) 

Motivated by the discussions in in Refs. flT4],|2"8|| , we choose the following "standard" value 



of the vector coupling constant^ 

G v = 0.5 Gs = 10.12 GeV -2 . (8) 

It should be stressed that the vector and axial-vector terms cannot simply be ignored in 
the effective Lagrangian ([J), as it is often done. These terms are necessary for correctly 
describing the vector meson properties P7[^5| , for adjusting the nucleon axial charge g& p0| , 



etc. In context of the present study, the vector interaction is important, because it generates 
a net repulsive contribution in asymmetric matter, i.e. when the numbers of quarks and 
antiquarks are not equal. 



1 This value is somewhat different from the one used in our previous paper [25], although the ratio 
Gv/Gs is the same. 
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Let us consider homogeneous, thermally (but not, in general, chemically) equilibrated 
quark-ant iquark matter at temperature T . Let ap t \ (bp t \) and a~p X (bp X ) be the destruc- 
tion and creation operators of a quark (an antiquark) in the state p, A , where p is the 
3-momentum and A is the discrete quantum number denoting spin and flavour (color indices 
are suppressed). By using the plane wave decomposition of quark spinors in Eq. (g), it can 
be shown |25| that quark and antiquark phase-space occupation numbers coincide with the 



Fermi-Dirac distribution functions: 



< a p,x a p,x >= n Pf = 


exp 


( Ep f - 
V T 




+ 1 


< b p,x b p,x,f >=^Pf = 


exp 


(E pf - 
\ T 




+ 1 



(9) 
(10) 



where Epf 
antiquarks: 



m. 



p 2 and pLRf, /i R f denote the reduced chemical potentials of quarks and 



PRf — Pf ~ GyPvf j 
~P>Rf = ~Pf + GyPvf ■ 



(11) 

(12) 



In our calculations we consider the chemical potentials \ij and ~fij as independent variables. 
The assumption of chemical equilibrium with respect to creation and annihilation of qq pairs 
would lead to the conditions 



The explicit expression for the vector density can be written as 

Pvf =pf-p f , 



(13) 



(14) 



where 



Pf 



v 



d 3 p 



• n 



(27T)3 P f " J (27T) 

are, respectively, the number densities of quarks and antiquarks of flavour / and v = 2N C = 6 
is the spin-color degeneracy factor. The net baryon density is obviously defined as 
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d 3 p _ 
^ n Pf 



(15) 



6 f 

The physical vacuum (pf = ~pf — 0) corresponds to the limit ripf = npf = . 

Within the NJL model the energy density and pressure of matter as well as the quark 
condensates psf contain divergent terms originating from the negative energy levels of the 
Dirac sea. As noted above, these terms are regularized by introducing the 3-momentum 
cutoff 9 (A — \p\), where 6 (x) = ~ (1 + sgnx). Then the scalar density is expressed as 

d 3 p rrif 



Psf = v 



n P f + n pf -6(A-p) . (17) 



(2tt) 3 E pf 

The energy density and pressure are obtained in a standard way from the energy- 
momentum tensor corresponding to the Lagrangian (fj). They can be decomposed into 
several parts as 

e = e K + e D + e s + e v + e FM + e , (18) 

P = P K + P D + P S + P V + Pfm + Po , (19) 

These expressions include: 
the "kinetic" terms 

e K = » E / 7^3 E Pf ( H Pf + n Pf) ' ( 20 ) 



3f J (2Tt) 3 E p/ 

the "Dirac sea" terms 



e D = -P D = -uJ2 J 7^ E P f G (A - P), (22) 



/ 

the scalar interaction terms 



(2k) 



es = -Ps = ^ L Y,Plf> (23) 
z / 



the vector interaction terms 



S 



ev = Pv = ^fT,Pv f (24) 
z / 



and the flavour mixing terms 



e FM = -Pfm = -AK J] p Sf . (25) 

/ 



A constant e = — Pq is introduced in Eqs. flUf ) and (|TUp in order to set the energy density 
and pressure of the physical vacuum equal to zero. This constant can be expressed through 
the vacuum values of constituent masses, m™ , and quark condensates, p$f. These values are 
obtained by selfconsistently solving the gap equations (§) in vacuum, i.e. at npf = npf = . 

For a system with independent chemical potentials for quarks (jif) and antiquarks (p*) 
one can use the thermodynamic identity 

e = J2(HPf+~PfPf)-P + * T - (26) 
/ 

Then one can obtain the standard expression for the entropy density, 

J^y [ n Pf lnn Pf + (1 - n pf ) In (1 - n pf ) + n pf -> n pf \ . (27) 
By using Eqs. <^), @-(^) one can also show that the differential relation 



dP = Y, [Pf dpf + p f dp f ) +sdT (28) 
/ 

holds for any thermally (but not necessarily chemically) equilibrated process. 



III. SYMMETRIC QUARK ANTIQUARK MATTER 

In this section we study the multi-quark-antiquark systems at zero temperature. Let us 
consider a symmetric system with equal numbers of quarks and antiquarks for each flavour. 
This requirement enforces the chemical potentials of quarks and antiquarks to be equal, 

~P f = Pf- (29) 

In this case the net vector density is automatically zero for each flavour, pvf = 0. The net 
baryon number and electric charge are also zero. Such systems are especially interesting 
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because the contribution of repulsive vector interaction, Eq. ([24]) , vanishes in this case. 
One can view such systems as a compressed meson gas where mesons are melted to their 
elementary constituents, quarks and antiquarks. The single particle states of quarks and 
antiquarks in such a system are schematically shown in Fig. 1(b). To understand this picture 
one should simply realize that antiquarks are holes in the Dirac sea. Then one can imagine 
that a certain number of quarks from the negative energy states are collectively excited 
into the positive energy states, producing an equal number of holes. Due to the presence 
of valence quarks and antiquarks the mass gap will be reduced. In this situation one can 
expect the appearance of bound states. 

To characterize the flavour composition we introduce the strangeness fraction parameter 

= N s + N- S 

Ts N u + Nz + N d + Nz+N s + N s ' 1 ' 

where iV/(/) is the number of quarks (antiquarks) of flavour / . For simplicity we consider 
only the isospin-symmetric mixtures where N u = and N n = iVj. 
Fig. 2 shows the energy per particle, e = e/p 

tot) clS " 

function of total density of valence 

quarks and antiquarks, p tot = (Pf + Pf) ; f° r different r s . At low densities e tends to the 

/ 

sum of the constituent quark and antiquark masses in vacuum, weighted according to r s , 

e(ptot - 0, r„) = m™(r s ) = (1 - r„) m™ + r s m™. (31) 

With growing density, e first decreases due to the attractive scalar interaction. At higher 
densities e starts to increase, approaching slowly the limit of ideal ultrarelativistic Fermi gas. 
For each r s one can see the appearance of a nontrivial minimum corresponding to a bound 
multiparticle state with e m j n (r s ) < m™(r s ). The density of qq pairs in these bound states 
varies between 0.7 and 1.4 fm~ 3 , depending on r s . It is easy to see [25[] that the minimum 
in e at any r s corresponds to zero pressure. Thus, finite droplets of such matter could be in 
mechanical equilibrium with vacuum. 

A more detailed behaviour of e in the plane p u +d — Ps is shown in Fig. 3(a) (here 
Pu+d = Pu + Pd = 2 p u ). The condition of fixed r s corresponds to a straight line with slope 
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r s / (1 — r s ) starting from the origin. By inspecting the figure one can notice the valley of 
local minimaP] starting from e = 0.482 GeV per particle in the pure ss system (r s = 1) and 
descending to e = 0.304 GeV per particle in pure uu + dd matter (r s = 0). There is no 
potential barrier on the way from r s = 1 to r s = 0. Thus, a droplet with any r s ^ will 
eventually "roll down" in the state with r s = 0. Possible decay modes of the considered 
bound states are discussed in Sect. VII. 

The constituent masses of u and s quarks as functions of ptot are shown in Fig. 4 for 
different r s . As expected, the constituent masses decrease with density, signaling the gradual 
restoration of chiral symmetry. It is interesting to note that the constituent mass is more 
sensitive to the density of the same flavour /. For instance, in pure uu + dd matter, i.e. 
at r s = 0, the s-quark mass does not change much with density. In this case m s varies 
entirely due to the flavour-mixing interaction. At densities corresponding to the bound 
states (indicated by dots on the respective curves) the constituent masses drop significantly 
as compared to their vacuum values: m u ^ ~ 0.1m™ , m s ~ 0.3 m™. For zero bare masses, 
Too/ = 0, the constituent quark masses in bound states would be practically zero. This 
means that the bound states correspond to the Wigner phase where chiral symmetry is 
restored. 

Fig. 5 shows the chiral condensate (uu) which, in our notation (see Eq. @), coincides 
with the scalar density ps u (notice the minus sign on the vertical axis). Its behaviour is 
strongly correlated with the u-quark constituent mass (compare with Fig. 4). As before, 
the variation of the (uu) condensate is mainly sensitive to the density of valence w-quarks and 
antiquarks. At r s = 1, when this density is zero, the change in the uu condensate is caused 
by the flavour-mixing interaction. The behaviour of the (ss) condensate is qualitatively 
similar, but due to the larger bare mass, too s — 140 MeV, its density dependence is weaker. 

The u and s chemical potentials are shown in Fig. 6. It is interesting that initially they 



2 This valley is shown by the dotted line in Fig. 3(a). 
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drop and then rise with density. One can easily understand this behaviour by analyzing the 
explicit expression for pf, 

fif = Jmj+p$ f , (32) 

1 /3 

where ppf = (Qir 2 pf/iy) is the Fermi momentum of quarks with flavour / (antiquarks 
have the same Fermi momentum in symmetric matter). At low densities, when ppj <C uif 
and fif ~ rrif, all chemical potentials decrease with density together with the respective 
constituent mass rrif (see Fig. 4). At high densities, when ppj 3> rrif, the chemical potential 

1 /3 

grows with density as in the free relativistic Fermi-gas, /if ~ ppf oc pj . In fact, this 
nontrivial behaviour of /i/ is responsible for the first order phase transition discussed below. 

IV. STRANGE QUARK MATTER 

Let us turn now to quark matter with nonzero net baryon density at T = 0. We assume 
that only valence quarks are present, i.e. the density of valence antiquarks is zero for each 
flavour (~pj = 0). In other words, this means that all levels in the Dirac sea are filled up, and 
additionally some levels in the Fermi sea are occupied by the valence quarks. This situation 
is schematically shown in Fig. 1(c). In contrast to the symmetric qq matter discussed above, 
now the symmetry between positive and negative energy states is lost due to the presence 
of repulsive vector interaction. 

Fig. 7 shows the energy per baryon e as a function of baryon density, ps = kj2Pf- 

f 

Different curves correspond to different r s , which in this case is simply the ratio of the 
strange quark density to the total density of all quarks. At ps — » the energy per quark 
tends to the corresponding vacuum mass given by Eq. fl31~|). With growing density both the 
attractive scalar and repulsive vector interactions contribute to e. It is interesting that at 
r s < 0.7 the attractive interaction is strong enough to produce a nontrivial local minimum 
at a finite ps- In the pure u, d matter (r s = 0) this minimum is unbound by about 20 MeV 
as compared to the vacuum masses of u and d quarks. On the other hand, it is located at a 
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baryon density of about 1.8 po, which is surprisingly close to the saturation density of normal 
nuclear matter. Of course, the location of this minimum depends on the model parameters. 
Nevertheless, one can speculate that nucleon-like 3-quark correlations, not considered in 
the mean-field approach, will turn this state into the correct nuclear ground state. 

When r s grows from to about 0.4, the local minimum is getting more pronounced 
and the corresponding baryon density increases to about 3.2 po- At larger r s , the minimum 
again becomes more shallow and disappears completely at r s ~ 0.7. At 0.2 < r s < 0.6 the 
minima correspond to the true bound states, i.e. the energy per quark is lower than the 
respective vacuum mass. But in all cases these bound states are rather shallow: even the 
most strongly bound state at r s ~ 0.4 is bound only by about 5 MeV per quark or 15 MeV 
per baryon. Nevertheless, the appearance of local minima signifies the possibility for finite 
droplets to be in mechanical equilibrium with the vacuum at P = 0. It is natural to identify 
such droplets with strangelets, which are hypothetical objects made of light and strange 
quarks [pl|-|36|. Similar multiparticle bound states made of nucleons and hyperons were also 
discussed p?|-p9|. 

It should be emphasized here that /3-equilibrium is not required in the present approach 
(see the discussion below). That is why our most bound strangelets are predicted to be more 
rich in strange quarks (r s > 1/3) than in the approaches assuming /^-equilibrium |3ll , |32| , |34 



which give r s < 1/3. As a result, these strangelets will be negatively charged^. Indeed, the 
ratio of the charge Q to the baryon number B is expressed through r s as 

%=*e»-l£±-l£L = l (l_3r s ). (33) 
B 3p B 3p B 3p B 2 K s) K } 



For r s ~ 0.4 this gives Q/B ~ —0.1. In light of recent discussions (see e.g. Ref. |p6[ ) 
concerning possible dangerous scenarios of the negatively-charged strangelet production at 
RHIC, we should emphasize that the strangelets predicted here are not absolutely boundP], 



3 Negatively-charged strangelets have been also considered in Refs. |32.35]. 
4 The analogous conclusion has been made in Ref. [24]. 
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i.e. their energy per baryon is higher than that for the normal nuclear matter. Hence, the 
spontaneous conversion of normal nuclear matter to strange quark matter is energetically 
not possible. 

The behaviour of the energy per quark for arbitrary mixtures of light and strange quarks 
is shown in Fig. 3(b). One can clearly see that the valley corresponding to the local minima 
has always a positive slope in r s . There is no potential barrier separating the states with 
r s ^ and r s = 0. Thus, a strangelet with any r s ^ will freely roll down to the normal 
nuclear matter state with r s = 0. It is interesting to note that a barrier in r s may appear if 
attractive scalar interaction is arbitrarily enhanced in the strange sector ||35|| . 

Fig. 8 shows the constituent masses of u and s quarks as functions of baryon density. The 
dropping masses manifest again a clear tendency to the restoration of chiral symmetry at 
high densities. As before, the dots indicate the masses at the local minima in the respective 
energies per baryon shown in Fig. 7. These mass values are somewhat larger than in the 
symmetric qq matter discussed above (see Fig. 4). Note that the stronger is reduction of 
constituent masses the deeper are the corresponding bound states. For the metastable state 
at r s = 0, which is a candidate for the nuclear ground state, the masses of u and s quarks 
are equal to 0.3 and 0.9 of their vacuum values respectively. For the most bound state at 
r s ~ 0.4 the respective mass ratios are reduced to 0.15 and 0.6. 

The behaviour of the u and s chemical potentials is shown in Fig. 9. One can notice 
the differences as compared with symmetric qq matter (Fig. 6). In particular, the curves for 
different r s do not intersect. Due to the additional contribution of the vector interaction, 
which is linear in ps (see Eq. ( |Tl~D where = ^Jm J + p F j at T = 0), the minima in \i u 
and \i a are less pronounced. At r s > 0.7 the curves have no minima at all. Therefore the 
chiral phase transition will be not as strong in this case as in symmetric matter. 
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V. SYSTEM ATICS OF BOUND STATES 



The properties of the multiparticle bound states discussed in preceding sections are 
summarized in Figs. 10-11. Fig. 10 shows the binding energy per quark, m™ c (r s ) — e m i n (r s ), 
where m™ c is the energy at ptot = and e m j n (r s ) is the energy at a local minimum, both 
taken for a given r s value. To avoid misunderstanding, Figs. 11(a) and 11(b) show e m j n (r s ) 
separately. Although the absolute minimum of e m i n corresponds to r s = 0, the maximum 
bindings are realized at nonzero r s . For symmetric qq systems the maximum value of about 
90 MeV is reached for r s ~ 0.7. This is indeed very strange and very bound matter! For 
asymmetric systems, where p/ = 0, the maximum binding energy is much smaller, about 
5 MeV per quark at r s ~ 0.4. One should bear in mind, however, that in this case e m j n results 
from a strong cancellation between the attractive scalar and repulsive vector interactions. 
Therefore, it is very sensitive to their relative strengths. The results presented above are 
obtained for Gy = 0.5 Gs- For comparison in Figs. 10 and 11(b) we also present the model 
predictions for Gy = 0. In this case the maximum binding energy increases to about 30 MeV 
per quark and the corresponding r s value shifts to about 0.6. It is interesting to note that for 
Gy = the bound state appears even in the pure u, d matter. The corresponding binding 
energy is about 7 MeV per quark, i.e. 21 MeV per baryon. 

Figures 10 and 11 reveal some differences compared to the previous calculations in 
Ref. [2S|. This is mainly due to a different set of model parameters used there. In par- 
ticular, the present set of parameters gives higher quark constituent masses in the vac- 
uum, m™^ ~ 368 MeV and mj ac ~ 550 MeV, as compared to m™^ = 300 MeV and 



m vac _ MeV in Ref. p5| . This leads to an overall upward shift in the energy per 
particle by 30-60 MeV. On the other hand, inclusion of the flavour-mixing interaction low- 
ers the energy. As a result, all binding energies increase slightly. A local minimum in e 
appears now even in the pure u, d matter. 

The dots in Fig. 11 indicate the positions of some conventional mesons (a) and 
baryons (b). Their empirical masses are rescaled (by factor 1/2 for mesons and 1/3 for 
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baryons) and shown at r s values corresponding to their flavour compositions. By inspecting 
the figure, one can make a few interesting observations. According to Fig. 11(a), mesoballs 
lie lower in energy than conventional vector mesons, but higher than the pseudoscalar mesons 
(see discussion in Sect. VIII). As one can see in Fig. 11(b), conventional baryons are more 
bound than strangelets even at Gy = . This is an indication that baryon-like 3-quark 
correlations might be indeed very important in the baryon-rich quark matter. 

VI. QUARK MATTER AT FINITE TEMPERATURES 

In this section we study properties of deconfined matter at finite temperatures. The 
calculations for this case can be done by using general formulas of Sect. II with the quark 
and antiquark occupation numbers given by Eqs. (j9l)-(|10|) . 

First we discuss thermal properties of meso-matter where the chemical potentials obey 
the conditions (p9|). A typical behaviour of the equation of state for this case is illustrated 
in Fig. 12(a). It shows the pressure isotherms for symmetric qq matter with strangeness 
content r s = 1/3. One can clearly see a strong first order phase transition which is signaled 
by the appearance of isotherms with negative slopes, d p P < (spinodal instability). The 
corresponding critical temperature is about 100 MeV. Another important feature is that 
the zero-pressure states persist up to temperatures as high as 70 MeV. This means that a 
finite droplet of this matter which has cooled down to this temperature and still remaining 
at high density of qq pairs (around 4p i n this case), will be trapped in a bound state. At 
later times it will further cool down by emitting hadrons from the surface (see below). 

The critical temperatures for a phase transition and for the appearance of a bound state 
are shown in Fig. 13(a) as functions of r s . One can conclude that this dependence is rather 
weak: the first critical temperature changes between 90 and 110 MeV while the second one 
varies around 70 MeV. 

Now let us consider chemically-equilibrated quark matter at finite temperatures. In this 
case the chemical potentials of quarks obey the conditions (|T3|). Fig. 12(b) represents the 
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pressure isotherms for the case /x s = . This condition implies equal numbers of strange 
quarks and antiquarks, i.e. zero net strangeness. It is appropriate for fast processes where 
strangeness is conserved, e.g. in relativistic nuclear collisions. As before one can see a region 
of spinodal instability, d p P < 0, which is characteristic for a first order phase transition. 
However, this phase transition is much weaker than in the case of symmetric qq matter 
(note the different scales in Figs. 12(a) and 12(b)). The corresponding critical temperature 
is about 35 MeV in this case. The zero-pressure states exist only at temperatures below 
15 MeV. 

Generally, the equation of state of the chemically equilibrated quark matter is charac- 
terized by two quantities: net baryon charge, B, and net strangeness, S. Therefore, it is 
interesting to study thermal properties of this matter at S ^ . Such states can be reached 
in neutron stars. They can also be realized via the distillation mechanism accompanying a 



QCD phase transition in heavy-ion collisions [j40|. Leaving a detailed study of this question 



for a future publication, here we only present results which complement our discussion con- 
cerning Fig. 13(a). Fig. 13(b) shows the critical temperatures for the phase transition and 
for the bound states in the equilibrated matter as functions of r~ = S/3B One can see 
that both temperatures first grow with r~ and then drop to zero at r~ ~ 0.8. The maximal 
values of respectively 50 MeV and 30 MeV are realized at some intermediate r~ around 0.4. 
As demonstrated earlier in Fig. 10, this value of r~ corresponds to the most bound state of 
strange quark matter at T = 0. So we see an obvious correlation: the deeper is a bound 
state at T = the stronger is a phase transition at finite temperatures. 

It should be emphasized here again that the thermal properties of asymmetric baryon- 
rich quark matter are very sensitive to the relative strength of scalar and vector interactions. 
To illustrate this point we have calculated the phase diagrams for several values of Gy. 



5 We introduce this new notation in order to distinguish this quantity from r s defined in Eq. 
In chemically equilibrated matter at T = (no antiquarks) these two quantities coincide. 
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The results for p s = are presented in Fig. 14. It shows the boundaries of two-phase 



coexistence regions calculated by using standard Gibbs conditions po| . If we take Gy = 0, 
as in most calculations in the literature, the coexistence region becomes wider and the 
corresponding critical temperature increases to about 70 MeV. On the other hand, if one 
takes Gy = 0.65 G$, the phase transition becomes weaker: the critical temperature moves 
down to 20 MeV and the zero-pressure states disappear completely. The calculation shows 
that there is no phase transition at Gy > 0.71 G,s ■ It is interesting to note that in all cases 
this first order phase transition occupies the region of densities around the normal nuclear 
density po . For instance, at Gy = 0.5 Gs the coexistence region at T = extends from 
0.2 to 1.7 p - 

These results demonstrate that this chiral phase transition is rather similar to a liquid- 
gas phase transition in normal nuclear matter. The critical temperature and baryon density 
in the present case (T c ~ 30 MeV, ps c ~ Po) are not so far from the values predicted 
by the conventional nuclear models |^T| (T c ~ 20 MeV, pb c ~ 0.5 po) • One may expect 
that in a more realistic approach, taking into account nucleonic correlations, the "chiral" 
transition may turn into the ordinary "liquid-gas" phase transition. If this would be the 
case, one should be doubtful about the possibility of any other QCD phase transition of the 
liquid-gas type at a higher baryon density. At least only one phase transition of this type 
is predicted within the NJL model. 



VII. DISCUSSION OF DECAY MODES 

Let us discuss briefly possible decay channels of the novel states of quark matter described 
above. Naively one could think that the symmetric qq matter would be extremely unstable 
with respect to the annihilation of quarks and antiquarks. But, in fact, many annihilation 
channels are closed or do not exist in dense qq matter. From the first sight, one may think 
that two-pion annihilation, q + q — ► 7r + 7T, should be the strongest decay channel in the 
bulk (one-pion annihilation is not allowed by the energy-momentum conservation). But 
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simple arguments show that this might be not true. Indeed, the bound states of qq matter 
appear at such high densities when chiral symmetry is practically restored. In this case 
the pion looses its special nature as a Goldstone boson and inside this matter it should be 



as heavy as other mesons. Moreover, due to the Mott transition ||42|| , it is unlikely that 
conventional mesonic states survive in this dense medium. Therefore, one can expect that 
the hadronic channels of the qq annihilation simply do not exist in the bulk. 

Another possible decay channel is the strong flavour conversion, ss — > uu or dd. In 
principle, this process can proceed through the one-gluon (color octet) intermediate state. 
However, it may be shown that due to color neutrality the corresponding matrix elements 
vanish after summation over the color indices. Within the present approach, the 4-fermion 
interaction terms are diagonal in the flavour space (see Eq. (0)) and, therefore, cannot change 
the flavour. Strong flavour conversion is possible only through the 6-fermion flavour-mixing 
interaction, ss — > uu + dd . This channel is only open when fi s > \x u + fid- By inspecting 
Fig. 6 one can see that this condition is never fulfilled. Therefore we conclude that strong 
flavour conversion is not possible or at least strongly suppressed in the bulk. 

In a finite droplet, hadronic annihilation, cff + qp — > hi + h 2 + is certainly possible 
at the surface, if fif + fif > + rrih 2 + ... . It is clear from Fig. 11(a) that the emission 
of pions in annihilation of light qq pairs is energetically more favorable than the emission 
of heavier mesons. Moreover, many annihilation channels of strange quarks and antiquarks 
are simply closed e.g. s + s — > K + K . As a result of the pion emission, the strangeness 
content of the daughter droplet will increase^} Finally, it will contain predominantly the ss 
pairs, similarly to a system composed of 0-mesons ("0-ball"). 



Further hadronization may be slowed down by several reasons ||25|| . According to Fig. 10, 
at r s ~ 1 the energy available in the annihilation of a ss pair is fi s + jx s ~ 0.95 GeV. This 
is lower than the threshold energies of the hadronic states 0(1020) and 7^^(990). The only 



6 Similar processes leading to strangeness distillation have been considered in Refs. [ 33 , 40 1 . 
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open channels in this case are p7r(910) and 37r(420). Their partial width in the 0-meson 
decay is about 0.7 MeV. If this would be the dominant decay mode, one could expect a life 
time for these 0-balls of about 280 fm/c, i.e much longer than the typical duration of a 
heavy-ion collision. 

In strange quark matter (without antiquarks), flavour conversion is only possible 
through weak decays. As follows from Fig. 9, at densities corresponding to zero 
pressure, the condition /i s > ji u holds. This means that weak processes of the types 
s — > u + e~ +V e , s + u — > u + d are allowed. Since there is no local barrier at any r s (see 
Fig. 11(b)), a system produced initially at some r s ^ will subsequently reduce this r s value 
by a cascade of weak decays. As a result, the system will roll down along the line of local 
minima shown in Fig. 11(b) (see also Fig. 3(b)). Finally, all s quarks will be converted into 
light u, d quarks. Two steps of this conversion process are shown schematically in Figs. 15(a) 
and 15(b). 

This picture is very different as compared to the one based on the MIT bag model. 
Within that model the pressure and energy density are given by simple expressions 

Pmyt = Hd ~ B, e M iT = e id + B , (34) 

where B is a bag constant, P;d and e;d are the pressure and energy density for a mix- 
ture of ideal gases of u, d and s quarks with constant (bare) masses, m u> d = m 0u d ~ 0, 
m s = m os — 150 MeV. By properly choosing B one can always get a zero pressure point 
and, accordingly, a minimum in the energy per baryon. Because the quark masses are 
kept constant, the condition ji s > m s will be first satisfied at a relatively low baryon den- 
sity oc rriQ S . At higher densities a certain fraction of s quarks will always be present in a 
/3-equilibrated matter. 

In contrast, in the NJL model the s quark mass is a function of both baryon density and 
strangeness content. As one can see from Fig. 9, at any given r s the condition fj, u ^ = fi s can 
be satisfied only at sufficiently high baryon densities which correspond to positive pressure. 
On the other hand, at the points of zero pressure we always have /i Ut d < fJ> s - Therefore, weak 
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decays will proceed until a strangelet reaches r s = (see Fig. 15(c)). 



VIII. CONCLUSIONS 

In the present paper we have investigated properties of deconfined matter at different 
densities of quarks and antiquarks, not necessarily constrained by the conditions of chemical 
equilibrium. All calculations are carried out within the SU(3)-flavour NJL model including 
scalar, vector and flavour-mixing interactions. We have demonstrated the possibility of 
strongly bound states in symmetric qq systems consisting of equal numbers of quarks and 
antiquarks ( "mesoballs" ) . The maximal binding energies of mesoballs, of about 90 MeV per 
particle or 180 MeV per qq pair, are realized for flavour compositions with about 70% of 
ss pairs. These systems remain bound up to the temperatures T < 70 MeV. The lifetimes of 
mesoballs may be long enough due to the suppression of annihilation into hadrons in the bulk. 
The model predicts a strong first order phase transition in chemically nonequilibrated meso- 
matter, i.e. at zero net baryon density. The critical temperature of this phase transition is 
in the range of 90-110 MeV depending of the relative abundance of ss pairs. As discussed 



in Ref. [25], formation of mesoballs in high-energy heavy-ion collisions may be observed 
through the event-by-event analysis of tt, K, spectra. In particular, narrow bumps in the 
hadron rapidity distributions may be generated by hadronizing mesoballs. 

By using the same model, we have also investigated the equation of state of chemically 
equilibrated deconfined matter at various temperatures, baryon densities and strangeness 
contents. The model predicts the existence of loosely bound, negatively-charged strangelets 
with maximal binding energies of about 20 MeV per baryon at r s ~ 0.4 . Similarly to 



Ref. [Q, no absolutely stable strange quark matter has been found. It is shown that prop- 



erties of baryon-rich quark matter are very sensitive to the relative magnitude of the vector 
and scalar interactions. At the standard value of vector and scalar couplings, Gy/Gs = 0.5, 
the metastable bound states of chemically equilibrated matter exist at T < 15 MeV, while 
at Gy = this temperature increases to 40 MeV. 
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The first order chiral phase transition in equilibrated baryon-rich quark matter is much 
weaker as compared to the symmetric qq matter. In the case of zero net strangeness the 
critical temperature is in the range of 30 MeV and the critical baryon density is around p . 
We believe that this phase transition is reminiscent of the ordinary liquid-gas phase transi- 
tion in nuclear matter. One can use the present model to study the chiral phase transition 
at nonzero net strangeness, which is appropriate e.g. for neutron star matter. We have 
found that the maximal critical temperature T c ~ 50 MeV is reached for the ratio of net 
strangeness to baryon charge S/B ~ 1.2 . 

It should be emphasized that all calculations in this paper have been made with constant 
values of coupling constants G$, Gy, K . We think that, this is a reasonable approximation 
at intermediate densities and temperatures. On the other hand, the asymptotic freedom 
of QCD requires all effective interactions to vanish at high densities and temperatures. 
Therefore, in a more realistic approach the coupling constants should decrease with density 
and temperature. Hopefully, this behaviour can be simulated by introducing finite range 
form factors. 
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FIGURE CAPTIONS 

FIG. 1. Schematic picture of energy levels of quarks and antiquarks in vacuum (a), symmet- 
ric qq droplet (b) and pure quark droplet (c). Occupied quark states are shown by 
full dots while open dots represent holes (antiquarks). R denotes the radius of a 
droplet, m and V are, respectively, constituent mass and vector potential of quarks. 
The boundaries of vacuum mass gap are shown by dashed lines. 

FIG. 2. Energy per particle e in symmetric qq matter vs. total particle density at different 
values of strangeness fraction r s (shown near the corresponding curves). Points 
indicate local minima of e. po — 0.17fm~ 3 is normal nuclear density. 

FIG. 3. Contours of energy per particle e (in GeV) for symmetric (a) and asymmetric (b) 

matter in the p u+ d — p s plane. Dotted lines represent local extrema of e . Shading 

shows regions with negative pressure. 
FIG. 4. Constituent masses of u and s quarks in symmetric qq matter as functions of total 

particle density. Figures in the box show values of strangeness fraction r s . Dots 

correspond to minima of energy per particle at given r s . 
FIG. 5. The same as in Fig. 4, but for condensate densities of u and s quarks. 

FIG. 6. The same as in Fig. 4, but for chemical potentials of u and s quarks. 
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FIG. 7. The same as in Fig. 2, but for asymmetric u, d, s matter without antiquarks. 
FIG. 8. The same as in Fig. 4, but for asymmetric quark matter. 
FIG. 9. The same as in Fig. 6, but for asymmetric quark matter. 

FIG. 10. Binding energies per particle in symmetric qq matter (solid line) and asymmetric 
quark matter as functions of strangeness fraction r s . Dotted line shows the results 
of calculations when the vector interaction is switched off (Gy = 0). 

FIG. 11. Minimal energies per particle in symmetric (a) and asymmetric (b) matter as func- 
tions of strangeness fraction r s . Dotted lines show energy per particle in the limit 
of zero particle densities, Eq. (131]) . Different parts of the solid line in the lower 
panel correspond to metastable (AB and CD) or bound (BC) states. Triangles and 
squares show masses of lightest mesons and baryons devided by the total number of 
constituents (two in mesons and three in baryons). Symbol < K > represents the 
spin averaged mass of K and i^*(892) mesons, i.e. (3m^. + m^)/8. Dashed line in 
the lower plot shows the results in the limit Gy — > 0. 

FIG. 12. Pressure isotherms for symmetric qq matter (a) and equilibrated quark matter with 
fx s = (b). Temperatures are given in MeV near the corresponding curves. Bound- 
aries of spinodal regions are shown by the dashed lines. 

FIG. 13. Critical temperatures for existence of bound states (P < 0) and phase transition 
(d p P < 0) in symmetric qq matter (a) and equilibrium quark matter (b) as functions 
of strangeness fraction. 

FIG. 14. Phase diagrams of equilibrium quark matter with zero net strangeness at various 
ratios of vector and scalar coupling constants. 
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FIG. 15. Schematic pictures of energy levels (shown by shading) occupied by light and 
strange quarks in cold quark matter at different strangeness fractions r s . Hatched 
boxes represent the constituent quark masses. Arrows show weak decay processes 
s — > u + e + V e and s + u — > u + d. 
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